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\^ ' based on the detection of broad Ha emission. The sample covers the range of black hole mass 



Abstract. We present the emission properties of a sample of 3 579 type 1 AGN, selected 



10^ < Mbh/M© < 10^-^ and luminosity in Eddington units 10"^ < L/LEdd < 1. Our main 

^ ' results are: 1. The distribution of the Ha FWHM values is independent of luminosity. 2. The 

I \ observed mean optical-UV SED is well matched by a fixed shape SED of luminous quasars, 

O . which scales linearly with broad Ha luminosity, and a host galaxy contribution. 3. The host 

galaxy r-band (fibre) luminosity function follows well the luminosity function of inactive non- 

^ \ emission line galaxies (NEG), consistent with a fixed fraction of ^ 3% of NEG hosting an 

AGN, regardless of the host luminosity. 4. The optical-UV SED of the more luminous AGN 
shows a small dispersion, consistent with dust reddening of a blue SED, as expected for thermal 
thin accretion disc emission. 5. There is a rather tight relation of vL^(2 keV) and broad Ha 

^ ■ luminosity, which provides a useful probe for unobscured (true) type 2 AGN. 

o 
o 
in 

^ ■ 1. Introduction 

p^ . We examine the emission properties of a new sample of low z broad line AGN, aimed to be an 

extension of the SDSS quasar catalog (QCV, Schneider et al. 2010) to low luminosities. This 
contribution summarizes the main results obtained on the distribution of the broad line FWHM 



(§3), the host galaxies (§4), and the AGN spectral energy distribution (§§5-7). Full details of 
the methods used, and a discussion of the implication of the results, are given in the main article 



C^ . (Stern & Laor, accepted). 

2. Sample Creation 

The type 1 (Tl) sample is selected from the SDSS 7*^ data release (DR7, Abazajian et al. 2009), 
based on the detection of a broad Ha emission line. The parent sample includes 232 837 spectra 
with S/N > 10 and 0.005 < z < 0.31. The broad Ha emission is identified by modeling the 
nearby narrow emission lines, stellar absorption features and featureless continuum. We find 
broad Ha emission in 3 579 objects (1.5%), which constitute the Tl sample. 

The optical spectra in the Tl sample are supplemented by nearly complete near IR 
photometry from 2MASS (Skrutskie et al. 2006, detection rate 97%) and UV photometry from 
GALEX (Martin et al. 2005, 93%). We also add X-ray photometry from the ROSAT (Voges 
et al. 1999) survey, available for 43% of the Tl objects. The measured broad Ha luminosities 
(LbHa) and FWHMs (Av) and the photometric luminosities are available electronically. 
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Figure 1. The distributions of the Tl (dark gray) and QCIV (light gray) samples in the broad 
Ha luminosity vs. FWHM plain. Crossing lines mark decades of the implied Mbh and L/LEdd- 



3. How broad do the broad lines get? 

Various theoretical models of the broad line region (BLR) and the AGN ionizing continuum 
predict an upper limit on Av, beyond which the BLR does not exist (Nicastro 2000, Elitzur and 
Shlosman 2006, Laor & Davis 2011). The models differ in the dependence of the maximal Av 
on luminosity. In Figure 1, we plot the LbHa vs. Av distributions of the Tl objects and the 
8185 QCIV objects where the broad H/3 is available (measurements from Shen et al. 2008). 

The Tl sample extends down by two orders of magnitude in LbHa compared to QCIV. The 
Eddington ratio L/Lsdd (= 10^ ^ind black hole mass Mbh (= lO^M©) are derived from LbHa 
and Av, in a manner similar to Greene & Ho (2005). It can be seen that the Eddington limit 
(/ = 0) sets a minimum Av with increasing LbHa ioi LbHa > 10^^*^ erg s"-*^. This, together with 
the rarity of AGN with m > 9.5, leads to a decrease in the range of observed Av values with 
increasing LbHa- A similar convergence of the range of Av for the Mg II and C IV lines at the 
highest continuum luminosities was noted by Fine et al. (2008, 2010). 

There is a steep decline in the number of objects with Av > 10 000 km s~^, at all luminosities. 
This decline is seen clearly in Figure 2, which plots the Av distribution of the combined Tl + 
QCIV sample at different Lboi- The distributions are remarkably similar, showing a roughly 
linear decline of log(dA^/dlog Av) vs. Av, or equivalently diV/dlogAv oc g-^vAvo^ with 
Avq ^ 2700 km s~^. The origin of this similarity is not clear, as the Av distributions should 
be set by the distribution of the m and / values. Either these distributions somehow lead to 
a Av distribution which is independent of LbHa? oi" it may imply that Av, rather than m and 
/, sets the observed Av distribution through some unknown mechanism. Note the effect of the 
Eddington limit at Av < 3 000 km s~^, which increases the minimal Av and the peak position 
with increasing Lboi- 



4. How does the AGN fraction change with galaxy luminosity? 

Kauffmann et al. (2003), and various follow up studies, found that AGN tend to reside in 
massive hosts. We therefore examine the Tl fraction as a function of host luminosity. 
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Figure 2. The FWHM distribution of the 
Tl+QCIV sample, at different Lboi bins. 

Figure 3. The spectral type distribution of 
galaxies in the parent sample, as a function of 
luminosity (r band). Non-Tl classifications 
are from Brinchmann et al. (2004). The 
'Tl host' line is corrected for the AGN 
contribution. The lower panel shows the 
number of objects per 0.25 decade in uLj^. 
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Figure 3 presents the classification oi z < 0.2 galaxies in the parent sample (§2) as a function 
of the spectroscopically measured uLjy{6166A). The type 2, Star Forming Galaxy (SFG), Non- 
Emission line Galaxy (NEG) and composite classifications (based on the equivalent width and 
BPT ratios of the narrow lines) are taken from Brinchmann et al (2004). At low luminosities 
>90% of the SDSS galaxies are SFG, while at high luminosities the NEG dominate. The 'Tl 
host' line is derived from the z/L^y(6166A) of the Tl objects, corrected for the AGN contribution 
(see §5 below). Remarkably, the Tl host fraction follows well the NEG fraction, or equivalently, 
their galaxy luminosity distribution is similar. This may suggest that Tls reside in NEG, as 
also implied by their similar concentration indices and bulge/total light ratios (Kauffmann et 
al. 2003). The fraction of NEG hosting broad line AGN at the level detectable in this study is 
~ 3%, independent of galaxy luminosity. 



5. How does the AGN spectral energy distribution (SED) depend on luminosity? 

Figure 4 compares LbHa with the observed 2.2 fim - 2 keV continuum luminosity, in the Tl 
sample. In the luminous bins in the FUV, NUV and 3940A panels, the slopes are consistent 
with a linear relation, in all other bands and luminosities the slopes are flatter. Since the FUV 
occurs close to the peak position of the AGN SED (e.g. Zheng et al. 1997), the linear relation 
between z/L^y(FUV) and LbHa suggests that LbHa provides a good estimator of Lboi- Also, the 
linear relation suggests that the mean covering factor of the BLR is independent of luminosity. 

In all bands excluding X-ray, there is a transition from a steeper slope to a flatter slope 
with decreasing LbHa- The transition LbHa increases with wavelength, up to a maximum of 
2^ 10^^ erg s~^ at 1.2/im. The transition occurs due to the host contribution, which peaks at 
1.2/im. 

The transition between AGN and host dominance of the SED can be seen in Figure 5. The 
left panel plots the mean SEDs in the LbHa bins, i.e. each line connects the mean values of 
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Figure 4. The LbHa vs. observed (AGN+host) z/L^^ from the X-ray to the NIR, in the Tl 
sample. The X-ray vLy are shifted upward by 1 decade for presentation purpose. The large 
black dots and associated error bars mark the mean and dispersion of vLy in half-decade wide 
LhYioL bins. The solid lines connect the means of adjacent bins, with the local slopes written 
underneath. Dashed lines depict a linear relation, normalized to the highest LbHa bin. 



different continuum wavelengths for the same LbHo; (Fig- 4). Note the transition from an AGN- 
dominated UV peaked SED to a galaxy-dominated NIR peaked SED with decreasing LbHa- As 
the net AGN SED, we take the Richards et al. (2006, hereafter R06) SED derived from luminous 
quasars in the SDSS. We scale the R06 SED to the UV in the highest LbHa bin, and rescale 
by the relative LbHa of the lower seven bins. The scaled R06 SEDs are the putative net AGN 
SEDs, assuming the mean SED remains fixed with luminosity and scales linearly with LbHa- 
In the middle panel of Fig. 5 we plot the residuals derived by subtracting the AGN from the 
observed SED, and in the right panel we plot the mean SEDs of inactive SDSS galaxies, with 
^-distributions matched to the ^-distributions of each LbHa bin. 

The residuals all appear to have a typical galaxy SED. This justifies (qualitatively) the simple 
scaling prescription we used for the net AGN SED. Note the significant excess luminosity in the 
Kg and H bands in the most luminous residuals, and the non-linear z/Ljy(X-ray)-LbHQ; relation 
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Figure 5. The luminosity dependence of the mean Tl objects SEDs at 2.2/im - 1528A. Left 
Solid lines mark the observed SEDs for the eight LbHa bins (mean LbHa noted). Dashed lines 
mark the net AGN SED, derived by scaling the R06 quasar SED by the mean LbHa of the bin 
with the same color. Middle Residuals of the subtraction of the the scaled R06 SEDs from 
the observed AGN SEDs. Right Mean SEDs of z-matched inactive galaxies (mean z noted). 
Dotted hues mark bands with < 70% detections. 



(Fig. 4), which imply that the linear scaling of the net AGN SED is applicable only to the 
optical-UV region. 

In order to quantitatively evaluate the AGN SED scaling law, we define the scaling index /3, 
where Lagn(A) oc i^^j^^. A /3 independent of wavelength A indicates the mean AGN SED shape 
remains fixed with luminosity. From the qualitative analysis shown in Fig. 5, we expect /3 ^ 1 for 
A in the optical or UV. In Figure 6, we plot the implied mean host luminosity Lhost = ^^z/ — ^AGN 
in the SDSS u and z bands, for different LbHa and /3. As the expected mean L^ost of ^ certain 
^bHa bin, we use the mean luminosity of type 2 AGN with the same z distribution as the LbHa 
bin. The type 2s are taken from the sample of Brinchmann et al. (2004). 

The top panel of Fig. 6 shows that in the u band, /3 < 0.9 is ruled out, as the implied residual 
is negative for LbHa > 10"^^ erg s~^. It is remarkable that the simplest scaling, ^5 = 1 in both 
bands, leads to the closest match of the mean type 1 host to the mean type 2 host. Together 
with the linear z/L^y(UV)-LbHQ; relationship at high luminosities (Fig. 4), this implies that the 
mean optical-UV net AGN SED has a fixed shape, and scales hnearly with LbHa- 
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Figure 6. A comparison of the implied mean 
Tl host luminosity for different LbHa bins, 
for different AGN SED scaling laws (black 
solid/dotted lines, /3 noted), with z-matched 
type 2 AGN. The comparison is made in the 
u band (3551 A), z band (8932A), and u — z 
color. The mean z of each LbHa bin is noted 
on top. Arrows mark a negative implied host. 



Figure 7. The Lx vs. LbHa relation, as a 
probe of true type 2 AGN. Dots mark Tl 
objects with FbHa > 10"-^^-^ erg s"-*^ cm~^. 
Red triangles mark mean upper limits in 0.5- 
dec LbHa bins. The best power law fit and 
dispersion are shown as black and gray lines 
(slope noted). Arrows mark the upper limits 
on LbHa of five true type 2 candidates. 



6. Is the lack of BLR emission in true type 2 candidates significant? 

True type 2 candidates are AGN in which the BLR is not observed, despite apparently 
unobscured X-ray emission. In Figure 7 we compare the upper limits on LbHa in five true type 2 
candidates from Shi et al. (2010) and Tran et al. (2011), with the Tl sample Lx(= i/L^(2keV)) 
vs. LbHa relation. Upper limits on LbHo; are derived from the expected Av (based on the 
published Mbh and L/LEdd? §3) and the flux density near Ha. The absence of a broad Ho: in 
NGC 4594, NGC 4698, and NGC 3147 is not significant. In the other two objects the expected 
^bHo; is well above the upper limits, and these two objects appear to be true type 2 AGN. 



7. What drives the dispersion in the optical- UV SED of AGN? 

In §5 we have shown that the mean SED of broad line AGN is well reproduced by the sum of 
the mean SED of luminous quasars, scaled down by LbHa? and a host contribution. Here, we 
examine the dispersion in the shapes of individual SEDs. 

In Figure 8, we show the distributions of z/L^y(FUV)/LbHQ; and optical slope aopt as a function 
of ^bHa- The z/L^y(FUV)/LbHa distributions for logLbHa > 40.5 have a very similar shape. All 
distributions have a narrow peak, with an extended tail towards low z^L^y(FUV)/LbHQ; values. 
The narrow peak suggests there is a small dispersion in the covering factor of the BLR. In the 
logLbHo; = 40.5 bin, host contribution to the UV broadens the distribution and increases the 
mean value. 
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Figure 8. The Tl objects z/L^(FUV)/LbHa 
and aopt distributions as a function of LbHa- 
Solid (dashed) lines represent bins of width 1 
(0.5) dec in LbHa- The mean LbHa are noted. 
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Figure 9. The o^opt vs. z/L^(UV)/LbHa of 
the Tl objects in the log LbHa — 43.5 bin. 
Also shown are the best fit (dashed line) and 
expected relations for different extinction laws 
(solid lines, marks every 0.1 in £'b-v)- 



The aopt distributions are redder at lower luminosities, due to the host emission. The 
log LbHa — 43.5 bin is least affected by the host contribution, and has a narrow peak and 
extended red tail, similar to the z/L^y(FUV)/LbHck; distribution, and to the continuum slope 
distribution found by Richards et al. (2003) on luminous quasars. There is therefore a small 
dispersion in z/L^y(FUV)/LbHQ; and o^opt in the net AGN emission, based on the objects least 
affected by the host contribution. The similar shape may indicate a common origin of the 
dispersion, such as dust reddening. 

The dust reddening scenario is explored in Figures 9 and 10. In Fig. 9 we plot o^opt vs. 
z/L^y(FUV)/LbHa foi" the logLbHa = 43.5 bin. The two independent ratios are highly correlated, 
as expected if dust extinction contributes to the dispersion of the AGN SED. An SMC/LMC 
extinction law (using the Pei 1992 formulation) is preferred over a Milky Way or Gaskell & 
Benker (2007) extinction law, as found by Hopkins et al. (2004) on luminous quasars. In Fig. 
10, we plot the 2.2 /im - 2 keV mean SED for different o^opt bins at the highest LbHa bins. In 
the logLbHo; = 43.5 panel, the difference between SEDs increases from the NIR to the UV and 
disappears in the X-ray, also consistent with dust extinction. A similar behavior is observed in 
the log LbHa — 43 panel, together with a weak anti-correlation of NIR luminosity and o^opt- 
The latter trend originates from the range of host luminosities, where a relatively luminous host 
increases z/L^y(1.2/im) and decreases o^opt- 

The reddening origin for the dispersion in the net AGN SED implied by Figs. 9 and 10, 
already at o^opt > — 1, indicates that the intrinsic AGN SED dispersion is even smaller than seen 
in Fig. 8. A rather uniform and blue intrinsic optical-UV SED is expected if AGN are powered 
by a thin accretion disc, which emits locally close to a blackbody. 



8. Conclusions 

We present and analyze a new sample (Tl) of 3 579 broad Ha selected AGN from the SDSS 
DR7, with log LbHa = 40 - 44, which spans m = 6 - 9 and / = -3 - 0. We add UV (GALEX), 
IR (2MASS), and X-ray (ROSAT) luminosities to form the mean SED. The main results are: 

(i) The Ho; FWHM velocity distribution dN/d log Av is independent of luminosity and falls 
exponentially with Av. The origin of this distribution remains to be understood. 

(ii) The observed mean 9000A-1500A SED, as a function of LbHa, is consistent with a sum 
of the mean SED of luminous quasars, which scales linearly with LbHa, and a host galaxy 
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Figure 10. The mean IR to X-ray SEDs and optical spectra as a function of aopt, foi" Tl 
objects with LbHa — 10^^'^ and 10^^ erg s~^. The value of Lx is shifted upward (+1 dec) for 
presentation. The error bars denote the errors in the mean Lx and z^L^y(FUV). 



contribution. 

(iii) The host galaxy r-band luminosity function of Tl objects follows the NEG luminosity 
function, with a relative normalization of ^ 3%, suggesting that the host of broad line 
AGN are NEG, and the AGN probability of occurrence is independent of the host mass. 

(iv) The dispersion in the optical-UV SED in luminous AGN (logLbHo; ^ 43), is consistent with 
reddening. This indicates the intrinsic SED of AGN is blue, with a small dispersion, as 
predicted from thermal thin accretion disc models. 

(v) The LbHa versus Lx correlation provides a useful probe for unobscured narrow line AGN. 
It can be used to test if the absence of a broad Ha, in X-ray detected AGN, is significant. 
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